Analysis of the Phase- Controlled Loop 
with a Sawtooth Comparator 

By A. JAY GOLDSTEIN 

(Manuscript received October 18, 1961) 

Because of the recent interest in phase-controlled oscillators, a discussion 
of the phase-controlled loop with a sawtooth comparator is presented. The 
main emphasis is on finding the pull-in range of the loop. A companion 
paper in this issue ( Ref. 4 ) deals with applications and shows how design 
parameters can be obtained from results developed here. 

I. INTRODUCTION 

The phase-controlled oscillator has evoked much interest in recent 
years. Some of its applications are to synchronism in television, 1 ' 2 syn- 
chronization to a harmonic of a crystal oscillator, 3 elimination of jitter 
in pulse code modulation, 4 tracking filters, etc. 

The general phase-controlled oscillator loop is given in Fig. 1. The 
incoming signal and the variable oscillator have the same free-running 
frequency o> c . The phase comparator has as its output some function / 
of the phase difference <p e = <p t — <fo . As examples of f{<p e ) we have 

the linear case : f(<p*) = <p e 

the sinusoidal case: f(<p.-) = sin ip e 

the sawtoothed case f(<Pe) = <?<• for — - < <p„ < - 

(see Fig. 2): f(<p e + nd) = /(*.*) for n - • • • -1,0,1, 

The output of the phase comparator passes through a filter whose im- 
pulse response is h(t). The output of the filter v(t) controls the variable 
oscillator according to the equation 

(103 
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Fig. 1 — The general phase-controlled loop. 
Thus, the frequency of the controlled oscillator is 

«. + *? = <*e + CCV(t). 

at 

In a companion paper in this issue, C. J. Byrne 4 discusses the engineer- 
ing origins and applications of the sawtoothed comparator and shows 
how design parameters can be obtained from the results of this article. 

This article is primarily concerned with finding the pull-in range of 
the loop. This is defined precisely in Section III. Briefly it is the maximum 
asymptotic (in time) value of the mistiming dtpt/dt for which the slave 
oscillator eventually synchronizes or locks to the in-put frequency. All of the 
literature cited in the references deals with this problem for the case of 
a sinusoidal or linear phase comparator. The linear case is easily solved 
since the resulting differential equation is linear. (See in particular 
Labin 5 for a detailed discussion. ) In the sinusoidal case the differential 
equation of the system is nonlinear. Only in the cases of no filter and an 
ideal integrator has the equation, up to the present, been solved in closed 
form. See Labin'' for an excellent discussion of the no-filter case. In order 
to handle the nontrivial filter, many authors have used methods of 
phase plane analysis. 6 ' 7,8 Phase plane analysis is restricted to the prob- 
lem of capture range in which the mistiming and phase error are zero 
for negative time, and the mistiming is constant for positive time. This 



T>e) 




Fig. 2 — The sawtoothed phase comparator characteristic. The phase error 
<f> c is difference between the input and output phases of the loop. 
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kind of analysis gives only upper and lower bounds for the capture 
range and is restricted to a lag filter (Fig. 3). For an RC filter (R-> = 0), 
Barnard 8 shows how phase plane analysis can give exact results. 

To obviate the mathematical complexities, people have resorted to 
making various hypotheses about the nature of the solution of the non- 
linear differential equation. These assumptions are based upon physical 
intuition and gross behavior observed in the laboratory. Different as- 
sumptions have led to different approximate solutions for the capture 
range. Moreover, they deal primarily with the lag filter, since it leads 
to a second-order differential equation while a more general filter gives 
a higher-order differential equation. 

The loop equation when expressed as an integral equation is 

^f = -a f f[<p e (t'))h(t - t') eft' + ^ - avo(t). 
at ->o at 

It is surprisingly tractable for the sawtooth comparator, and the pull-in 
range can be computed for any filter. Fig. 4 shows the excellent agree- 
ment between theory and experiment for the lag filter. These experi- 
mental results were obtained by C. J. Byrne. 

To obtain our results, we too must make an assumption. While the 
assumptions other authors have made deal with the behavior (in steady 
state) when far outside the pull-in range, ours deals with the behavior 
just outside of the pull-in range (see Section 4.4). This hypothesis is 
easily verified experimentally and has been so verified by C. J. Byrne 
for a representative selection of RC filters. 

A brief description of each section follows. 

Section II gives the basic integro-differential equation of the loop. 

Section III defines the lock and pull-in range. The former is called by 
some the pull-out range. The lock range is the maximum frequency 
difference that the loop can lock to. It is given by 

w/. = cxf nmx _H(0) 

r 



o — vw 

Rl 



<*R2 

k 



1 



i"he integral compensating or lag filter. The normalized time con- 
stants are t, = a (ft, + It t )C and t 2 = aR»C. For an RC filter t 2 = 0. a = (V. F. O. 



Fig. 3 - I 

...ants are n . 

output frequency shift )/(V. F. 0. input voltage) 



606 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1962 




40 60 100 



Fig. 4 — The relative pull-in range. For critical damping (r* + l) 2 /4 = n and 
for the RC hlter ll 2 = 0. 

where H(0) is the dc gain of the filter and /max is d/2 for the sawtooth 
comparator. 

Section 4.1 gives the solution of the basic loop equation. This solution 
is the sum of ( 1 ) the solution of the linear phase comparator problem, 
(2) a series of step functions, and (3) a series of damped exponentials. 
The solution is obtained by representing the phase comparator function 
as the sum of the phase difference [giving (1)] and a series of translated 
unit step functions [giving (2) and (3)]. 

Section 4.2 gives the steady-state solution when not captured. In this 
case the output of the phase comparator is a periodic function whose 
period for a fixed filter depends on the asymptotic relative mistiming 
(Fig. 5). By examining this non-capture situation we obtain the pull-in 
range. We observe that in non-capture state the period and relative 
mistiming must correspond to a point on a curve typified in Fig. 5. 
Hence a relative mistiming lying below the minimum point of the curve 
corresponds to a capture or synchronized situation, and the height of 
the minimum gives the ratio of pull-in to lock range (the relative pull-in 

7p). 

Section V gives all the explicit design formulae for the lag filter. For 

the special case of the RC filter (R 2 = in Fig. 3) an explicit formula 

for relative pull-in can be given, namely 



7p 
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tanh £ (aRiC - £)~* (olfr 



C>07 



(aRiC < i 



1). 



In all other cases we must find the roots of a transcendental equation by 
numerical approximation methods. 

Byrne 4 gives graphs of the results of Section V for the lag filter. These 
are graphs of relative pull-in (Fig. 13), noise bandwidth (small signal) 
(Fig. 7), figure of merit (relative pull-in/noise bandwidth) (Fig. 15), 
and maximum loop gain (small signal) (Fig. 8). 

The noise bandwidth is a measure of the ability of the loop to reject 
small phase noise. More explicitly, the noise bandwidth N of a network 
is defined to be the bandwidth of that ideal low-pass filter which passes 
the same white noise power as the given network. 

There are many possible ways of defining a single measure of the 
performance of the system, depending on the particular application in 
mind. We have chosen the figure of merit y P /N, i.e., a large figure of merit 
implies high noise rejection and large relative pull-in. 




T = PERIOD 



Fig. 5 — Relative mistiming u m /u L in a non-synchronized steady state vs the 
period T of the comparator output, (a) no filter, (a) and (b) overdamped loop and 
(c) underdamped loop. 
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For small phase deviations of the input, the comparator can be con- 
sidered linear. We can then discuss the gain of the loop as a function of 
the frequency of the phase deviation. The maximum of the loop gain is 
denoted by Y. In some applications Y is restricted by stability considera- 
tions to be less than unity. 

Section VI is devoted to the derivation of several interesting asymp- 
totic results for the lag filter. A simple formula is obtained for the 
relative pull-in for large values of the filter time constants. It is also 
shown that if the maximum loop gain is allowed to have a fixed value 
greater than unity, then, by appropriate choice of the time constants, 
arbitrarily large values of the figure of merit can be obtained. 

This work could not have been completed without the aid of M. 
Karnaugh who suggested the problem, E. G. Kimme who proved that 
the sawtooth comparator is a continuous approximation to the original 
discrete sample data system, C. J. Byrne whose experimental work con- 
firmed the formulae derived here, D. E. Rowlinson who constructed the 
contour curves from the computer data, and R. D. Barnard with whom 
many fruitful discussions were held. 

II. THE BASIC LOOP EQUATION 

We obtain an integro-dirTerential equation for the loop by noting that 
the output of the filter can be written as a convolution plus initial condi- 
tions 

v(t) = I f[<P,(t')]h(t - dt' + v (t) 

where v (t) is the filter output due to residual charges and fluxes in the 
filter at time zero. v (t) damps out exponentially in all filters of interest. 
Substituting this into (1) and replacing <p by <pi — ip e we obtain 

^ = -a f f[<p e (t')Mt - dt' + *£ - av (t). (2) 

at Jo at 

In order that the derivations which follow not be unduly complicated 
by inessential parameters, we make the following normalizations 

x(t) = <p c {t)/f ma * (/ max = d/2) 

C(x(t)) =/(fp.(0)//-« 

<p(t) = <p,-(0//n.» • 
The normalized form of (2) becomes 

% = -a f C[x(t')]h(t - dt'+%- otvo(t)/f mx . (3) 

dt Jq dt 
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III. DEFINITIONS OF LOCK RANGE AND RELATIVE PULL-IN 

If the input frequency w c + dipi/dt is increased "very slowly" to a 
value which is not too large, the output frequency u> c + dtpjdt will 
follow it (i.e., be always equal to, or locked to, the input frequency). 
The maximum value of dipi/dt for which lock-in will occur is called the 
lock range and is denoted by co L . More precisely, u L will be determined 
from (1) when the maximum dc voltage v is obtained. This maximum 
value is clearly the product of a, / max the maximum value of the com- 
parator function /and //(0) the dc gain of the filter.* 

ul = af m „H(0). (4) 

Suppose that the input frequency is not increased slowly, but in some 
sudden or erratic manner. Suppose moreover that the input frequency 
approaches a limiting value, u>,„, the mistiming; i.e. 

i- dtpt 

lim —r- = 0) m . 

i -»oo dt 

In general, even if < <a m < u L (that is, we are in the lock range), the 
output frequency will not asymptotically lock to the input frequency 
(that is, be captured), but will be a modulated frequency. We define 
the relative ■pull-in range y p to be that normalized maximum frequency 
difference such that 

— ypfi>L < Hm -£ = u m < y p w L (5) 

{•Me at 



implies 



lim § = co,,, . (6) 

<-.» dt 



Notice that we make no restriction on how d-p/dl approaches to,,, , as 
long as | o)„, | < y p <a L . 

IV. DERIVATION OF RESULTSf 

4.1 Basic Equation 
Let 

< to < k • ' ' < tn < • • ■ 



* We shall use capital letters to denote the Laplace transform of the function 
denoted by corresponding lower-case letters. 

t From here on \vc are dealing with the sawtooth comparator. 



u(t) = 
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be all the instants (called discontinuity points) at which the phase 
difference x(t) crosses the discontinuity of C, i.e., 

limx(t n - A) = x{t n -) = 1 +2n' 

A-»0 
A>0 

where the first equality is a definition of x{t n -) and where n' is an 
integer dependent on n. Let 

a n = 1 if x is increasing at t n 

a n = — 1 if x is decreasing at t n 

a n = if x is stationary at t n . 

Using the unit step function 

[0 for t < 

} 1 for t ^ 

we can express C(x(t)) in the analytically useful form 

£WL > = *® - n. - t "Ml - W (7) 

2 2 y=o 

where w is an integer so chosen that this equation holds at t = 0. 
We note here for future reference that 

n 

x(t n -) = n ±h + £«;• (8) 

Substituting (7) into the loop equation (8), we obtain 

ldx = _af l x{t , )h(t _ t>) dt , + ano r h{l!) d( f 
2dt 2 J J o 

+ a E aj [ u(t' - tj)h(t - t') dt' 
i—o Jo 

+ lJ--b(«)/V. 



(9) 



olving this by Laplace transform methods we obtain 

ai e- sii aH(s) 



Yf \ - s3> (s) - [yo(0) + «7o(s)]//max , Wo ocH(s) 
i AW ~~ 2(s + aH(s)) " «8 + aff(«) 



y=o S s + aH(s) ' 
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Letting 



R(s) " X ^7 \ < 10 > 

s -f- aH(s) 



we have 



i -*(.)- _^w_. (11) 

S + aH{s) 

Note that sR(s) is the transfer function between input phase <p, and 
comparator output phase y e for the linear comparator case. r(t) is then 
the phase response at the linear comparator output due to a step in 
input phase. Since applications will require the system to synchronize 
to a step in phase, we will assume that r(t) —> as t— » ». 

Using this equation and taking inverse transforms in the equation for 
X(s) we obtain 

w. = *M + no(1 + r(t)) + £ a . w( , _ ^ _ £ a . r( , _ ,. } (12) 

* ^ y=o y-o 

where 

v , , s*(s) - fo,(0) + «7o(«)]//«„ 
s + a//(s) 

.rz.(/) is the solution of the loop equation in the case of a linear com- 
parator function f(x) = x. 

Using the final value theorem 10 we have 



Xl( oo ) = lim.r L (0 = lim sX L {s) = lim <p'(t) 

(-»« B-*0 t-,00 



x L (<*>) = 



aH(0) 

2a>„,/d <a m 



(13) 



aH(O) u L 

From (12) and (7) we have for the comparator output 

C(x(t)) x L (t) A . 

5 — = —5 2^ajr(t - tj) + n r(t). (14) 

In a steady-state condition this reduces to 

C(x(t)) = ^-o^ajrit- tj) (15) 

Oil, — °° 

where the w KO term vanishes because of the remarks following the 
definition of R(s). 
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4.2 Steady-State Solution When Not Captured 

When we are not locked and in steady state, the output of the phase 
comparator will be a periodic function.* We give here a simplified 
heuristic derivation of the steady-state periodic solution. A rigorous 
derivation is easily obtained using the heuristics as a guide. In steady 
state, the normalized comparator output y(t) = C(x(t)) will be periodic 
with a period which we will call T. In a given period there may be many 
discontinuity points tj ; let us suppose there are k. Then assuming we 
are in steady state, we can write 

m _i_ t /«=••■ -1,0,1, ••• ,.,m 

where 

= To < Ty < ■ ■ ■ < 2Vi < T. 
These relations are illustrated below. 

\t„k \tnk+l | I' ' 'I |*nfc+fc-l \kn + D* 

nT + t nT +T t + T ■■■ nT + T k -i + r (n + l)T + r 
|< T >| 

The a„'s will be periodic in steady state and we let 

fn = ••■ -1,0,1, • • • , m 

It is no restriction to assume a time shift so that t = 0. Then, let 

/ = mT + m (0 < u ^ T) (18) 

and combine the above three equations with (14). We obtain 
y{t) = C[x(t)] = C[x(mT + u)] 

= o)„,/u) L — 2 2 a n k +i r(mT + u — <„*+,) 

*— 1 m 

= Wm /« L - 2 IV-1,- 5>[(m - n)T + u - T,]. 



(The second summation has the upper limit m because r(t) = for 
/ ^ 0.) Letting) = TO - n, we obtain 

y(t) =^-2Ei,E r(jT + u- T,-). (19) 



Oil, , = j = 



* A mathematical proof is not at hand. Indications of its truth arc given in 
Benes 9 and experimental observations confirm this. 
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Let us define a periodic function 

P (t,T) = iS r(t+JT) <°<'* T ) (20) 



\p(t - nT,T) (nT <t^ (n + l)T) 



or 



p(t,T) = £r(f +jT). 

(r(t)) = Ofor / ^ makes p(t,T) a well defined function.) With this 
definition, the normalized steady-stale comparator output, when not 
locked, can he written 

2/(0 =22-2 lUtfO- r„r). (2i) 

CO/, , = 

The expression for p(t,T) is familiar to those in the field of sample 
data systems.* Though superficially formidable, it can be expressed in 
closed form quite easily for the only important class of the filter transfer 
functions //(.s), namely rational functions. In that case R(s) is a rational 
function too. Hence r(t) is a linear combination of exponentials of the 
form t'"e ' (real part of negative). Then p(t,T) for < t £ T is a 
linear combination of geometric series, each of the form 

z(l) = E (/ + jT) m e P(t+jT) 

, = 

= -£i>*' +jr) (22) 

<tp'" j=0 

,m at 

d e 



ttpm i _ e ar 

This steady-state solution consists of a constant term 2a> m /daH(0) 1 
which is the normalized steady-state output for a linear phase com- 
parator plus a linear combination (with coefficients ±1) of time trans- 
lates of the function p(t,T), which is periodic of period T. The derivation 
shows that every steady-state periodic solution of the loop equation has 
the form of (21). 

Equation (21 ) hides several pitfalls. These are: 

1. We must have | //('/) | ^ 1. Hence only certain T and 7\ are 
admissible. 



* It is the response of a filter li(s) to an input ^ *(' + JT). 

i— * 
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2. Are the solutions represented by (21) physically realizable? 

3. Are the solutions represented by (21) stable with respect to small 
noise perturbations? 

These three topics are grouped under the title Boundary Conditions 
and will be discussed following a discussion of the pull-in range. 

4.3 Relative Pull-in 

From the definition of T, y(T- ) = ±1 and by an appropriate choice 
of r in (16) (if co m > 0) we may assume y{T) = 1. Then from (21) 



a— i 



- B = l + 2i; A iV (T - Ti ,T). (23) 

u L «'=o 

Now the minimum value of u m > for which we have a non-constant 
periodic steady-state stable solution is by definition y p u L , hence 

y p = 1 + 2 min £ A iV {T - Ti , T). (24) 

i=0 

where the minimum is taken over all T and over all steady-state solu- 
tions satisfying conditions 1, 2 and 3 above. 

4.4 Boundary Conditions 

4.4.1 Discontinuity Point Condition 

y(t), being the normalized phase comparator output, satisfies — 1 ^ 
y(t) ^ 1. Also y(t'— ) = ±1 if and only if for some n and i, t' = Ti + 
nT, or y(t) is stationary at t' (i.e., y'(t') = and y at t' is increasing if 
y(t') = -1 or decreasing if y(t') = 1). These are equivalent to 

Y,Mp(t' - T i9 T) - P(T - Ti,T)) = 

if and only if if = nT + T { or y(t') - y(T) is stationary at t'. This 
restriction will be called the discontinuity point condition. 

To analytically determine whether this condition is satisfied, in a 
general case, is clearly very difficult. For the case of the lag filter we 
can solve the problem analytically but must rely on an experimental 
fact. C. J. Byrne has found experimentally, in a large class of RC filters, 
that there is just one discontinuity per period T, i.e., the k in (21) is 
one. We will call this the Experimental Hypothesis. Thus 



y(t) = ^ - 2p(t,T) (25) 

0)L 



and 
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y p = 1 + 2minp(r,T). (26) 

T 

In the section on the lag filter we show that if 

ViT'X) = mm p(T,T) (27) 

T 

then p(l,T') satisfies the discontinuity point condition. Thus if p(t,T') 
is realizable (it is — see below) and is stable under noise (we do not 
know, but have some evidence — see below) then 

7p = 1 + 2p(T',T') (28) 

for the lag filter. 

4.4.2 Realizability Condition 

Does there exist, for each of the steady-state functions represented 
in (21 ) satisfying the discontinuity point condition, a corresponding 
input function ip(t)1 That is, are the y(t) in (21) physically realizable? 

In Appendix A we prove realizability for any filter but not in quite 
the form stated above. We do the following: 

(a) A particular input (p(t) = 2<a m /d is injected. 

(b) The loop is broken at the output of the phase comparator. 

(c) Into the filter, at this point, is injected a voltage which asymptot- 
ically has the form (21). 

(d) One shows that the output of the phase comparator has asymp- 
totically the same form. 

(e) In steady state the loop is closed. 

4.4.3 Non-Synchronous Stability 

Are the solutions stable? By this we mean: Will a steady-state solution 
be thrown into synchronism by a "small" noise? In formal terms, we 
suppose that a solution y(t) has a discontinuity point, say / shifted by 
noise to / + A . Each of the following discontinuity points h , t 2 , • • • , 
/„,••• is shifted to t x + Ai , t 2 + A 2 , • ■ • , L + A„ , • • • . It suffices for 
our purposes that the (/„ -f- A„)'s be asymptotically periodic (i.e., the 
noise sends us into another periodic solution and not into synchronism). 
The host we have been able to prove is that 



lim 



dA„ j 

rfAo U„-o J 



This has been done for the lag filter using the experimental assumption 
that k = 1 and that 7" — e < T < T', for e sufficiently small, where 
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T' is given in (28). Now it would suffice for stability to show that A„ is 
hounded for A sufficiently small, but the above does not imply this, for 
all it says is that 

A„ = cA + e„A 2 

and we do not know that e„ is bounded. 

V. LAG (INTEGRAL COMPENSATING) FILTER 

5.1 General Results 

This section gives all the explicit formulae for design procedures in 
the case of the lag filter (Fig. 3). We assume the experimental hypothesis 
(see Section 4.4) throughout this section. 

The transfer function of the filter is 

ff(.) = ks + L 



where 



Hence 



M+ 1 

h = (Ih + RjC. 
k = R*C 



R( \ = I = k8+ 1 

(S) s + aH(s) fis 2 + {ah+ 1)8 + a 

Pl + ^ 1 P2 + T 2 1 



Pl — P2 S — Pi Pi — P2 S — p 2 

where pi and p 2 are the roots of denominator of R(s). In particular, 
introducing the normalized dimensionless time constants 

n = ati , i = 1,2 



we have for the roots 



where 



Pi =i(a+ (-1)*'6) 

h 
a= (r 2 + l)/2 ^ i 



* The real or imaginary part of b is non-negative. 
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The denominator of B(s) can be written in the form 

s 2 + 2co„£s + co„" 
where 

to,, = ( a/i\ ) 
and £, the damping factor, is 

£ = ( T2 + l)/2 V7i = a/ (a 2 - I)' 1 )*. 
In this notation we obtain 

r(7) = 1 [-(a - 5 - 1) exp(-(a - 6)</<i) 
26 

+ (a + 6 - 1) exp(-(o + 6)J/<i)]. 

Because ;•(<) is a linear combination of exponentials, we can easily sum 
the infinite series for p(l,T), obtaining 

, , exp[-(o + 6)77'] 

+ (a + 6 - 1) f — 7 — T-rn 

1 - exp[-(a + «)i?]_ 

where 77' = J//j and y = T/t x are dimensionless time variables. 
To obtain y p using the results of (27) we must find 

min p( T,T) 

T 

or the roots of 

= <lp(v,v) 

Differentiating the expression for p( 77,77) we obtain 77 9* and 

sinh 2 (a — 6)77/2 _ (a — 6)(a — 6 — 1) ,„_v 

Hiiih 2 (a + 6)77/2 ~ (a 4- 6)(a 4-6-1) 

or 77 = x . And upon using the addition formula for the hyperbolic sine, 
we have 

tanh 077/2 6 tanh 677/2 _ 9 a 2 4- 6" — a ^ _ , ., . > 

I tanh 077/2 " 2a — 1 

£ tanh 677/2 

which defines c, or 77 = « . 
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Use of the quadratic formula gives 



tanh ari/2 

r tanh 677/2 




= c + Vc 2 - 6 2 = Ciiafi) (32) 



or 77 = 00 .* In special cases considered it was found that the minimum 
of p(rj,Ti) occurs at the first positive zero of its derivative (or at 77 = 00 ). 

5.2 Critical Damping 

From (31) we see that as b approaches zero (damping factor equals 
one), 

tanh 077/2 n a(a — 1) / nN ., . , 

-nr— = 2 -k r = c i( a >°), lf a > ] 

V/2 2a - 1 (33) 

77= oo, if I < a < 1 . 
Thus y p = 1 for 6 = and § < a < 1. 

5.3 ATo Filter and RC Filter 
The filter parameters satisfy 

g r 2 ^ r, 

which upon conversion to the a and b parameters become 

(a - l) 2 ^b 2 
and 



Equality holds in the first case, when Ri = or C = (i.e., there is no 
filter) and in the second case, when R 2 = 0, (i.e., a simple RC filter.) 

For no filter, a -\- b — 1 = or a — b — 1=0, and referring to (30) 
we have only 77 = °°. Thus min p(T,T) — and y p = 1. 

For the RC filter R 2 = 0, a = £, we obtain from (30) 77 5* and 
sinh brj = or 77 = 00 . If & is real, 77 = =c and 7 P = 1. If & is imaginary 

77 = mir/b m = 1,2, • • • 



* If the negative sign were used in the quadratic formula then 77 would be nega- 
tive (complex) when b was imaginary (real). 
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and we easily find that p{rmr/b,mir/b) is minimum at m = 1, giving 
finally 



tanh - 



u-r 



if T1 >1 



(34) 



if 



Tl ^ 



<i 



The results of these special cases are graphically summarized in Fig. 
6. (Also see Fig. 13 of Byrne, Ref. 4.) In the shaded area of Fig. 6 the 



ft 



0.5 - 



a. 





/a-b= i 






^a+b=i 








0.5 1.0 




(a) 



/ y 

(NO FILTER) Y ^ ^^-- 




^^ UNDERDAMPED 
S LOOP 




jf~«-P& 






(b) 



1.0 



Fig. 6 — In part (a) the parameters a and b are restricted to lie below and/or 
to the right of the polygonal curve. The heavy lines and the shaded area give 
values of a and b for which the relative pull-in is unity. In part (b) the same in- 
formation is given for the normalized time constants rj and r» . 
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relative pull-in is unity. This follows from the fact that the left-hand 
side of (31) is bounded below by 26* while 

2c - 2b = 2(a - b - l)(a - b)/(2a - 1) 
is negative in that region. Hence in (31) we must have 77 = <x> . 

5.4 Com-putatiorial Procedures 

Except in the special cases of no filter (Ri = 0) and the RC filter 
(R2 = 0), there is no simple way of computing the relative pull-in. We 
must solve (32) by an iterative procedure and substitute the result into 
the equation for p(rf,rj). If 77 is the solution of (32) or (33) we have a 
simpler equation for y p , namely 

y p = [1 — D seen 2 ai?/2]/tanh ar,/2 



where 



and 



D m (« - Do, - 6' 

Cl 2 — 6 2 



D = (a - \)/a (b = 0). 



An upper bound for r? is obtained from (32) and (33). Using the fact 
that tanh x < 1, we obtain 

. f2(tanh- 1 6/c 1 )/6 (6 * 0) ,„., 

' < \2/c 1 (6 = 0)" (35) 

A lower bound for r\ in the case b is real is obtained by using the in- 
equalities 

z — z /3 ^ tanh z ^ z. 

Using this in the equation for rj we have 

I , . v tanh br)/2) 
a V /2 - (ar,/2)'/Z ^ tanh a v /2 = } c ^ a > 0) ~ -5 I g ClV /2 

{ Cl (a,0) V /2 J 

giving the lower bound 

* The left-hand side of (31) is of the form b(x + l/x). For x positive this is 
bounded below by 26. 
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We note here for future reference that if b is imaginary, b = ib' then 
(35) implies the inequality 

^ 776' < 7T. (36) 

5.5 Discontinuity Point Condition for the Lag Filter 
To prove that this condition is satisfied, it suffices to show that 

d ± * for < t < T. (37) 

at 

For, since we may suppose 

</m = 1, 

it follows that if 

u(t') = 1 (0 < /' < 

then Rolle's theorem tells us that there exists a t" with t' < t" < T 
such that //(/") = 0. This contradicts (37). It suffices also to prove 
(37) for that T which minimizes p(T,T). 

Recall that we are assuming we have a lag filter and that /.• = 1 in 
(21) (experimental hypothesis). Assuming (37) false, we obtain from 
(29) after some calculation 



-[a+b)v'/2 , _ -(a+6)i}/2 



e -(a-MV/2 I _ g-(a-6)ij/2 

where rj minimizes p(u,u). Note that < t\ < t\. 
Case 1, b real. Then a > b and 

-(n+6).)'/2 

! = e - b "' 

e ~(a-b)r,' 72 



(38: 



> 


8 -H 




-(a+6)i)/2 




e -(a-6),/2 


> 


j _ -<a+6),/2 * 


1 _ e -(a-b) v l-l' 



Hence (38) is false. 



* If < x < y < 1, then j/f/ > x — 1/y — 1, for -x > — i/ implies j-;/ - j > 
xy — y; hence in factoring and dividing we obtain the desired inequality. 
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Case 2, b imaginary. Let b = ib', then (38) becomes 

-&V/2 + wtt = arg(e- (a+,V)V2 - 1). (39) 

Now the real part of e~ u * 2 — 1 is negative and the imaginary part 
is negative (since by (36), < b'rj/'Z < x/2). Hence the right-hand 
side is an angle in the third quadrant. But the left-hand side is an angle 
which can only be in the second or fourth quadrant, since 

< &V < b'v < v. 

Hence (39) is false, proving the discontinuity point condition. 

5.6 Small-Signal Properties of the Loop 

In this section we give formulae for design parameters of the loop 
when we are operating on the linear portion of the phase comparator. 
Then the closed loop transfer function Y is 

Y( s ) = *(** + *) 



ttti + (ah + 1)S + a' 

Restricting our attention to real frequencies and normalizing the fre- 
quency co by 

£2 = o)/a 

and recalling that 

we obtain 

T 2 2 fi 2 + 1 



F(fi)!" = 



1 (t s + l) 2 fl 2 + (1 - DPn) 1 ' 
With the phase shift 

(l + nntf)0 

B = — arctan 



(i - r,n 2 ) + r 2 (i + T 2 ) 2 n 2 

if denominator positive * 



\ir if denominator negative ' 

Important parameters for design are the maximum gain and the 
frequency and phase shift at which it occurs and the range of frequencies 
for which the gain exceeds one. Differentiating | Y(Q) | 2 and solving for 
its zero gives 

* The arctan is an angle in the first or fourth quadrant. 



PHASE-CONTROLLED LOOP 
2 



623 



(2n - l)/2n if t 2 = 0, n ^ \ 

{[1 + (2(r, - r 2 ) - Dr^/n 2 ] 5 - l}/r 2 2 if n - r a £ J 
[0 if r, - to ^ i 
Solving | 7(ft) | 2 ^ 1 gives 



ft 2 ^ ftj 2 



where 



(2(t, - t 2 ) - 1 . 



fti" = 







if r, - t 2 £ $ 

if Ti — T 2 ^ 5 



We also have the interesting inequality 

V2 fl ma x ^ fii 

with equality when t 2 = 0. The cases r 2 = and n — t 2 ^ \ are im- 
mediate. The ease n — t 2 ^ £ gives 

fi„,ax 2 = {[1 +n\k*t - MA/ 
_ fti 2 

= [1 + T 2 2 fii 2 ]* + 1 



proving the result in this case. 

We wish to emphasize that the maximum gain is unity if and only if 
n — t 2 ^ \. Peak gain = constant contours are given in Fig. 8 of Ref. 4. 

The 3 db point occurs at ft = ft} where 



r(fl0 I 2 = I 



from which we obtain 



where 



-2 N } 



ftr = B + (# + O 



£ = (to 2 + 2(r, - to) - 1)/2t 2 2 . 
The noise bandwidth N is defined by 



N 



-fl 



F(co) ! 2 r/co. 
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It can be evaluated in various ways, for example see Refi 10. One obtains 

.V = «*(1 + r 2 2 /ri)/2(r 2 + 1). 

In the no-filter case (t 2 = n = 0) and RC case ( T -> = R 2 = 0) we have 
N = ira/2. N = constant contours are given in lief. 4, Fig. 7. 

As discussed in the introduction, the figure of merit was chosen to be 
the ratio N/y p . N/y p = constant contours are given in Ref. 4, Fig. 15. 

VI. ASYMPTOTIC RESULTS 

In this section we obtain the asymptotic results stated in the introduc- 
tion. Since the derivations are tedious, the results are first summarized. 

From computer data, the contour curves of relative pull-in y p = 
constant with ordinate and abscissa the normalized time constants 

n = a(R x + R 2 )C 

t-> = aR%C 

seem to be asymptotic to straight lines for large values of the normalized 
parameters. (See Fig. 13 in Ref. 4.) This observation led to the con- 
jecture that for fixed y p and large t-> 

n = K{t 2 + 1). 

In Appendix we prove this and show that 

l/K = 1 - (1/7;, - y P )'\tanh^ y p f. 

With respect to the figure of merit (see Fig. 15 in Ref. 4), the following 
very important results are derived in Appendix B for the lag filter. Sup- 
pose the peak small-signal phase gain Y of the loop is restricted to be 
unity (it is always unity atdc). Then the maximum merit obtainable 
for filters giving the unity peak loop gain is 2.27. If, however, we permit 
a fixed peak gain greater than unity, we can have an arbitrarily large 
merit figure. This usually results in very poor transient response. More 
precisely, the following results are derived in Appendix B. Let us con- 
sider those lag filters for which the peak small-signal (phase) gain is 
fixed at F. Define M by 

M- = 1 - Y~ 2 

Then for a filter with normalized time constants r r and t 2 and normalized 
frequency ft = w/a, for which the loop has peak gain Y occurring at 
frequency £2 mn x , we have 
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and 

n = (il/V + 2r 2 + 1)/2(1 + M). 

Asymptotically for t 2 large we obtain for the noise bandwidth (with 
a = (r 2 + 0/2) 

2(1 - M)\ / ,_ , n/ _- 2N * 



AT/xa = ( 1 + " w - 1 / 4a + 0(a 

and for the relative pull-in range 

_. 2 y/M - 1 , „, -|v 

=^(^y + ° (( ^ ) - 

Thus the noise bandwidth decreases as a -1 while the relative pull-in 
decreases as a~\ Hence the figure of merit increase as a . 

The derivations of the preceding results are given in Appendices B 
and C. 



APPENDIX A 

Realizability of Steady-State Solutions 
Recall that (assuming d = 2) 

y& = -Wffu ~ 2 ^ *< £ r(t - Ti - nT) (40) 

where [see (13)] 



aH(0) ' 

Since we assume y(t) satisfies the discontinuity point condition 

fc— i 

= 1+2 
aH(0) 



= 1 +2T,A iP (T - Ti,T). 



Break the loop at the output of the phase comparator, inject y(t) 



* Two functions /(x) and g(x) satisfy /(a:) = 0(g(x)) if and only if | f(x)/g(x) | g 
constant < «= for x sufficiently large. 
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into the filter, and let the input phase be o) m t + x L — c (where c is 
defined below). The phase output of the oscillator is given by 



d<po(t] 



= a f y(t')h(t - t') dt' 



dt 'o 



and upon integrating once and substituting (40), 

- 2 £ Ai E ( f <xr(t" - T t - nT)h(t' - t") dt" dt'. 

,=0 n = ■'0 •'0 



By taking the Laplace transform of the double integral in the summa- 
tion and by using the relations in (10) and (11), we find 

*(0 =t^ (' f Ht")dt"dt' 

tl\\J) ■'0 J o 

- 2E^,(E u(t - Ti - nT) -£; r(t- T t - nT)\ . 

/=0 [n— n=0 J 

Now the remaining double integral is the integral of the step response 
of the filter and for large / is of the form H(0)t + c. Using this and the 
definition of y(t), we obtain for large t 

ft— 1 00 

Mt) ~ «J + c - 2 £ A t £ u(t - Ti - nT) - y{t) + x L . 

i=0 n=0 

Now using the discontinuity point condition and the representation of 
the comparator in (7) we find the comparator output is asymptotically 
y(t). Hence in steady state we may close the circuit without any dis- 
turbance. 

APPENDIX B 

Figure of Merit for Constant Peak Gain and Large Time Constants (Lag 
Filter) 

From Section 5.6 we have for the closed loop small-signal (phase) 
gain 

I Y(n) I 2 _ t 2 SI + 1 , v 

1 K n (t, + l) 2 fl 2 + (1 - fiVi) 2 * 
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Differentiating with respect to ft and equating the result to zero gives 

TtVtfnax + 2rX,„ - [2(7! - T 2 ) - 1] = 0. (42) 

We can also represent the square of peak gain Y~ as the ratio of the 
derivatives of the numerator and denominator of (41) evaluated at 



2 

y' 1 = r - 



(r 2 + l) 2 - 2n(l - tApc) 

2 

T2 



7 



., 2 - [2(r, - t 2 ) - 1] + 2ri«fl 



202 



iii:ix 



This, after using (42), gives 

F 2 = ^-r- . (43) 

Defining A/ ^ by 

m 2 = i - y- 2 , 

we have 

£ ilf < 1, since 1 ^ F < ». 
Also (43) gives 

nnLx = M. (44) 

Substitute (44) into (42) and solve for r, . Then 

r , . (MV + 2r, + 1)/2(1 - il/). 

Using this result in the formula for the noise bandwidth (Section o.fi), 
we have for Y constant and t 2 large (and hence a = (t 2 + l)/2 is 
large) 



.V 



We now turn to the problem of obtaining asymptotic expressions for 
the relative pull-in range for F fixed and greater than unity. 
We can rewrite the expression for t\ as 

ri = JML a > + 2(1 + M)a - *1±^ . (40) 



* If /(.r) = p(x)/q(x), then/'(x ) = implies /(x ) = p'(x )/q'(x ). One obtains 
this result by logarithmic differentiation of f(x). 
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Using the definition b 2 = a 2 — t x , (46) and the binomial expansion, 
we have for large a 

if M ^ | and 

b 2 = -3(a - i) (48) 

if M = 5. /n the following we suppose M 9^ \. Recall (31) that to find 
the relative pull-in we need the root of 

, /0 tanh 677/2 
tanh a-q/2 = c t - 



where 



and 



Hence 



Ci = c + (c 2 - 6 2 )' 
c = (a 2 + ft 2 - a)/(2a - 1). 
c = (2a 2 - a - n)/(2a - 1) 



c = 



= a[l — ti/o(2o — 1 )] = a — — (49) 



Also 



c 2 - b 2 = (a - IlY - (a 2 - 



T\ 



(50) 



giving 



Finally 



2 \ i T\ 



T2 



HfeH©'-©]- 
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c, = ( a - - 1 ] + (c- - b 2 ) 1 



(51) 
C -o 1 -=-+ 0(a" 2 ) . 

Setting 2 = aij/2, we have 

tanh 2 =» r tanh -2. (52) 

a 

We will show that for large a, 2 is small and then obtain an approximation 
to z by using a power series expansion for tanh z. First note that the 
derivative at zero of the right-hand side of (52) is 

*_ i-i +0(0 

a 2a 

which approaches 1 from below for large a. Also 

£ " + 0C6- 1 ) 

2A/ ' ' ' + 0(6- 



1 - M 



Hence | Ci/b | is bounded away from 1 (and greater than 1). 

A sketch of the curves of the two sides of (52) with the above two 
facts shows that 

lim 2 = 0. 

a— co 

Using power series expansions in (52) we obtain for large a 
3 b V \a/ 3 \a, 



or 



t = 3 ' - c ' /a 



1 — cib 2 /a 3 
1 



+ 0(0 
2AP~ 



-8-ii (53) 



1 - il/ 



+ 0(0 



= (3(1 / ))> a-'rO(a-). 
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From Section 5.4 the relative pull-in is 

y, = -L=° + D tanh (f) 
tanh grt 



where 



Then 



and 



since 



D= ((a- l) Cl - 6 2 )/(ci 2 -b 2 ). 
(d - a + 1) 



1 - D = 



ci — 6 2 /ci 



1 _ D = * ~ g + 1 (54) 

2c 



ci + 6 2 /ci = 2c. 
Using (49) and (51) we have 

1 - D - w + 0(a " !) - (55) 

We now obtain the asymptotic formula for 7 P by substitution into 
the formula for y p the approximations for D, 1 — D and the approxima- 
tion tanh (arj/2) ~ a-q/2 = z with z approximated as in (53). 



2(1 - 71/)' _j , n , _»s 

*" wr a +0(a,) 



= p((r 2 +l)/r 1 ) i + 0(r 2 -J). 



(56) 



APPENDIX C 

Relative Pull-in {J jag Filter, Large t\ and t 2 ) 
Assuming that for a large a 

n = 2Ka + L + 0(a _I ) (57) 

we obtain from the definition 

b 2 = a 2 - n 
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that 

Expanding the square root we obtain 



(58) 



and 



- = 1 - - + 0(0. (59) 

a a 



From (57) since 

a = (r 2 + l)/2 
we have 



t, r , , L + X 



12 = K + i^-L-2: + 0(r 2 - 2 ) (60) 

t 2 r 2 

From (50) we have 

'-'-©"-© 

and by using (60) we have 

^-»P-( g -- g )[, + ^ + gy i) -» I + o(o]. 

Using the binomial expansion 

(c 2 - b 2 ) 5 = (K 2 - ff)*[l +^ (L+ ^ ) i 2 ^" 1) ^+0(r 2 - 2 )]. (61) 

From (49) 

Tl 

c = a — — 

T2 

and using (60), we obtain 

c = a - K - *i±A + 0(r 2 " 2 ). (62) 

T2 

Then 
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c\ = c + (c 2 — b 2 ) 7 

= (a - K) + (K 2 - tf) J + 0(r 2 -'). 

Finally from (58) and (63) 

c, = d/(o - K) 
b b/(a - K) 

= i + <** ~ K)h + (a- 2 ) . 
a 



(63) 



(64) 



Letting z = ar}/2, (31) becomes 

tanh 2 = [1 + (/C 2 - K) h /a + 0(a~ 2 )] tanh ( 1 - K/a + 0(a~ 2 ) )z. (65) 

Using the addition formula for tanh (A + B)z and simplifying, we have 

tanh 2 * tanh (K/a + 0(a~ 2 ))z - [(K 2 - K)'/a + 0(a~ 2 )] tanh 2 

+ [1 + (K 2 - K) h /a + 0(a' 2 )] tanh (tf/a + 0(a~ 2 ))z = 0. 

We show that z/a approaches zero with a and use this to simplify (66). 
From (35) 

2 tanh -1 [1 + (K 2 - K) h /a + 0(<T 2 )] 



z < 



1 - K/a + 0(a- 2 ) 



i..l^#^„n + o<o 



1 - K/a + 0(a~ 2 ) 
Since 

lim In a/u = 

u-»°o 

we have 

lim z/a = 0. 



Returning to (66), we now have asymptotically 

~K 



. - iK — IV tanh 2 
tanh 2 4- I — 77 — —1=0. 



Solving for K we obtain 

\/K = 1 - — V- - tanh 2 2 2 . (67) 

Ltanh 2 J 
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Now the relative pull-in given in Section 5.4 is 

7„ = 5 — + D tanh z 

tanh z 

and we easily show that [using (54)] 

Ci — a + 1 



1 - D = 



2c 
[K 2 - K) 1 - K + 1 + 0(o" 



a - K + OCa- 1 ) 

= 0(0. 

Hence asymptotically for fixed z, 

7p = tanh z + 0(a _1 ). (68) 

Thus for given relative pull-in, the above gives us z and tanh z, and 
then (67) gives K from which (57) gives for large n 

r x = 2C(r 2 + 1). (69) 
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